Noise-induced synchronization for phase 
turbulence 

Hidetsugu Sakaguchi 
Department of Applied Science for Electronics and Materials, 
Interdisciplinary Graduate School of Engineering Sciences, 
Kyushu University, Kasuga, Fukuoka 816-8580, Japan 

February 8, 2008 

Abstract 

Phase turbulence is suppressed by applying common noise additively 
to the Kuramoto-Sivashinsky type equation, and the noise-induced phase 
synchronization is realized. The noise strength necessary for the suppres- 
sion of phase turbulence is evaluated theoretically. 

PACS: 05.45.+b, 05.40.+j, 02.50-r 

Recently, various noise effects to nonlinear systems have been studied. The 
response of a bistable system or an excitable system to a periodic force is en- 
hanced by the noise effect. The stochastic resonance improves signal detection 
by the superposed noise P3 Gl E]- Noise-enhanced entrainment among cou- 
pled oscillators is found experimentally in Belousov-Zhabotinsky reactions 0]. 
Frequency locking of noise-sustained oscillations is found in coupled excitable 
systems A small amount of noise may change a chaotic trajectory into a 

rather regular trajectory, and it is called noise-induced order [Jj. Common noise 
may induce complete synchronization for uncoupled chaotic oscillators [HHH]- It 
is called noise-induced synchronization. In this paper, we apply common noise 
for a modified equation of the Kuramoto-Sivashinsky equation. Without the 
common noise, the model equation exhibits phase turbulence. Noise-induced 
synchronization occurs and a spatially uniform state is observed owing to the 
common noise. 

The model equation is written as 

4> t = u> - r sin <j> - ^cj> xx - cj> xxxx + (j> 2 x + £(t), (1) 

where <j>(x, t) is a phase variable, u> > r is assumed, and £(t) represents spatially 
uniform Gaussian white noise satisfying 
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The model equation is equal to the original Kuramoto-Sivashinsky equation, 
when r = D = 0. The Kuramoto-Sivashinsky equation exhibits phase turbu- 
lence for fi, > 0, owing to the long- wavelength instability of the uniform state. 
For nonzero r, the phase motion of each oscillator is not uniform. This type 
of phase oscillator model has been used to study cooperative phenomema in 
coupled limit cycle oscillators Equation (1) represents a model of 

spatially coupled phase oscillators. Without the noise term, the spatially uni- 
form rotation is unstable and a phase turbulence is observed similarly to the 
original Kuramoto-Sivashinsky equation. 




Figure 1: (a) Stationary probability distribution P s (0o) for the Fokker-Planck 
equation at uj = 0.8, r — 0.5 and D = 0.7. (b) Average value of — rcos^o as a 
function of D. Below the dashed line (— r cos^o) = —0.0625, the uniform state 
is stable for uj — 0.8, r = 0.5 and [i = 0.5. 



Firstly, we study the stability of the spatially uniform state under the com- 
mon noise. The spatially uniform state (f>(x,t) — (f>o(t) satisfies the Langevin 
equation 

d<t>a/dt = uj — r sinc^o + (2) 
The corresponding Fokker-Planck equation is 

dP d d 2 P 

For D = 0, the stationary probability distribution is 

P s (M<x (4) 

uj — r sin <fio 

For nonzero D, the stationary probability distribution is expressed as [111 I12| 



P 3 (4>o) = exp{(uj^ -r+rcos(f) )/D}P s (0) 



J 2 " et-^-r^^^dip 
(5) 
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Figure 2: Snapshot profiles <f>(x,t) for Eq. (1) at co = 0.8, r — 0.5 and D = 
(a), 0.2 (b), 0.7 (c) and 2.5 (d). 



where P s (0) is determined from the normalization condition J w P s ((j) )d(j)Q = 1. 
This stationary probabiliy distribution approaches a uniform distribution P s = 
1/27T for D — > oo. The stationary distribution for u = 0.8, r — 0.5 and D = 0.7 
is displayed in Fig. 1(a). The peak of the probability distribution is located at 
<p ~ 7T/4 < 7r/2. 

The linearized equation of the Fourier amplitude 4>k with wave number k 
around the uniform state for Eq. (1) satisfies an equation 

d(j> k /dt = (-rcos0 o (t) + ^fc 2 - fc 4 )</> fe , (6) 

where 4>k = l/V~L (f>(x)e~ lkx dx, and 4>o(t) obeys Eq. (2). The long-time 
average of — rcos</>o(i) + /ifc 2 — fc 4 determines the stability for the perturbation 
with wave number k. For D = 0, the average value of cos^o is zero, since the 
distribution (4) is symmetric around 4> = Tt/2. For D — * oo, the average value 
of cos4>o(t) is also zero, since the distribution is uniform. For general parameter 
values of Z?, the average value of — rcos<po(t) is negative, since the probability 
distribution has a peak below n/2 as shown in Fig. 1(a), and cos</>o is positive 
for < 0o < 7r/2. We have calculated the average value of — r cos 4>o(t) using the 
stationary distribution (5) for various D's for uj — 0.8 and r = 0.5, and plotted 
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Figure 3: Average of Fourier amplitude Pw(k) = (\4>k\) at u — 0.8, r — 0.5 and 
D = (a), 0.2 (b), 0.7 (c) and 2.5 (d). 



the results in Fig. 1(b). The curve of (— r cos <f>o) has a characteristic form of the 
resonant one, although it is not easily expected from the distribution (5). The 
average value of — r cos </>o approaches for D — > and D — > oo and it takes a 
minimum value -0.08 at D ~ 0.7. The average value of — rcos^o(i) + A^ 2 — ^ 4 
takes a maximum — r(cos0o) + M 2 /4 at A; = ^/ /z/2, when fc is changed. The 
uniform state is linearly stable if — r(cos</>o) + /J. 2 /4 < 0. At cj = 0.8, r = 0.5 
and 13 = 0.7, the uniform state is stable if fj, is smaller than 0.4. For fi = 0.5, 
lik 2 — fc 4 takes a maximum value 0.0625 at k = 0.5, therefore, the uniform state 
is stable for 0.3 < D < 1.58. It implies that the common noise of intemediate 
strength can suppress long-wavelength fluctuations and stabilize the uniform 
state. Common noise of intermediate strength is most effective for the noise- 
induced synchronization, which is analogous to the stochastic resonance. 

We have performed numerical simulation to see the noise-induced synchro- 
nization for the modified Kuramoto-Sivashinsky equation. We have used the 
finite difference method with time step 0.0001 and space step 160/512 for a 
numerical simulation, and the periodic boundary conditions with system size 
160 is used. The parameter values are w — 0.8, r = 0.5 and fx = 0.5. The noise 
intensity D was changed as a control parameter. Figure 2 displays snapshot 
profiles of 4>{x) at D = 0,0.2,0.7 and 2.5. At D = 0, the spatio-temporal chaos 
is observed. At D = 0.7, the spatially uniform state appears. It corresponds 
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Figure 4: (a) Intermittent time sequence of a Fourier amplitude |0/s(i)| with 
wavenumebr k = 0.471 for ui = 0.8, r = 0.5 and D = 0.2. (b) Histograms of 
the probability distributions of \cf>k\ with wavenumber k — 0.471. The straight 
line denotes a power law with exponent -0.74. (c) Intermittent time sequence 
of a Fourier amplitude |c£fc(t)| with wavenumebr k = 0.942. (d) Histograms of 
the probability distributions of \4>k\ with wavenumber k — 0.942. The straight 
line denotes a power law with exponent -0.93. 



to the noise-induced synchronization. At D = 0.2 and D — 2.5, the spatial 
profile is not uniform, but fairly regular. The wavelength of the most dominant 
Fourier mode is 160/12, and the corresponding wavenumber is 0.471, which is 
close to a/ /i/2 = 0.5. The amplitude of spatial fluctuations changes intermit- 
tently in time. The intermittent time sequence of the Fourier amplitude 
with k = 0.471 for D = 0.2 is shown in Fig. 4(a). The amplitudes of the spatial 
fluctuations seem to differ very much for D = 0.2 and 2.5 in Figs. 2(b) and (d), 
but it is a problem of timing to take snapshots. The average amplitudes of the 
spatial fluctuations do not differ so much for the two parameter values of D, as 
shown in Figs. 3(b) and (d). 

Figure 3 displays the average of the Fourier amplitude: P w (k) — (\<fik\) for 
the four parameters of the noise intensity. For D = 0, the spectrum has a peak 
at k ~ vmT^ an d increases near k ~ 0. This is a characteristic spectrum of 
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phase turbulence. For D = 0.7, the uniform state is stable and the spectrum 
is almost zero. For D = 0.2 and D = 2.5, the spectrum has strong peaks at 
k ~ yj fj,/2 and its higher harmonics, but the long-wavelength fluctuations are 
suppressed in contrast to the case of D = 0. The suppression of the long- 
wavelength fluctuations makes the spatial profiles of <j>(x) regular, as shown in 
Fig. 2. 

The time evolution is intermittent, and <f)(x,t) is close to the uniform state 
<j>o(t), when the amplitude of spatial fluctuations becomes sufficiently small. 
Then, the Fourier amplitudes obey Eqs. (2) and (6) approximately. The variable 
4>o(t) is a random varibale, since Eq. (2) is a Langevin equation. Equation (6) 
represents a random multiplicative process, since the growth rate — r cos 4>o(t) + 
/ik 2 — k for the Fourier amplitude with wave number k is randomly fluctuat- 
ing. If the average value of the growth rate is positive (negative), the Fourier 
amplitude generally tends to grow from zero (decay to zero). Even if the av- 
erage value of the growth rate is positive, there are chances that the growth 
rate keeps negative in a certain time interval. The Fourier amplitude takes very 
small values for the time interval. Inversely, even if the average value of the 
growth rate is negative, there are chances that the growth rate keeps positive 
in a certain time interval. The Fourier amplitude takes non-small values for 
the time interval. This is a mechanism that a random multiplicative process 
induces an intermittent time sequence p3 EH EH ■ When the amplitude of spa- 
tial fluctuations becomes large, the nonlinear term acts and the growth of the 
spatial fluctuations is suppressed. Figures 4(a) and (c) display time sequences 
of Fourier amplitudes \(f> k \ with (a) k = 0.471 and (c) k = 0.942 at D = 0.2. 
The average growth rate is positive for wave number k — 0.471 and the average 
growth rate is negative for wave number k = 0.942. For both wave numbers, 
intermittent time sequences appear. Figures 4(b) and (d) display histograms of 
the probability distributions P(\4>k\) of \4>k \ for (a) k = 0.471 and (b) k = 0.942. 
The probability distributions obey approximately a power law for a rather large 
range of k, which is characteristic of multiplicative stochastic processes [141115) . 
The exponent of the power-law distribution is nearly -0.74 for k — 0.471 and 
-0.93 for k = 0.942. The exponents change with wave numbers. Similar type 
intermittent behaviors are also observed for D = 2.5. 

To summarize, we have found noise-induced synchronization for the modified 
Kuramoto-Sivashinsky equation. The phase turbulence is suppressed by the 
common noise. The mechanism for the noise-induced synchronization is clear in 
this model equation. That is, the stationary phase distribution of the uniform 
state is deformed by the noise and the uniform state is stabilized against the 
long-wavelength fluctuations. For r = as in the original Kuramoto-Sivashinsky 
equation, this type of noise-induced synchronization does not occur. The growth 
rate of the Fourier amplitude with the long wavenumber is decreased especially 
for an intermediate range of the noise sterngth. The average growth rate is 
positive for a finite interval of wave numbers, and the spatial profile becomes 
rather regular owing to the common noise. The fluctuation of the growth rate 
induces an intermittent time sequence and a power-law distribution for the 
Fourier amplitude has been observed. 
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